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Abstract. A new modeling way of describing the continuous absorption of 
electromagnetic (EM) radiation in dense partially ionized hydrogen plasma is tested in 
this work. It is shown that the obtained results give a possibility of calculating spectral 
^ 1 absorption coefficients which characterize the relevant absorption processes in partially 

■ ionized hydrogen plasmas with electron densities N e ~ 10 19 cm -3 and temperatures 

T f=a 2 ■ 10 4 -fT. A key parameter of the used procedure is determined empirically. 
The calculation method is applied to wavelength region 300nm < A < 500nm. The 
presented results can be of interest for dense laboratory plasmas as well as for partially 
ionized layers of different stellar atmospheres. 
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1. Introduction 

^ ■ By now, direct methods of determination of various plasma characteristics, based on 

quantum or classical statistical mechanics, have been developed only for practically fully 
ionized plasmas [THE]. In the case of dense partially ionized plasmas, where the density of 
neutral particles (atoms) is close to the density of positively charged particles (ions), such 
rigorous methods do not exist at present. Recently, in [9], this problem was discussed in 
connection with the transport properties of dense partially ionized plasmas. As for their 
optical characteristics, it is enough to remind that adequate calculation methods exist 
only for weakly and moderately ionized plasmas with electron density N e < 10 17 cm~ 3 . 
As it is well known, the influence of the neighborhood on an exited atom can be neglected 
in such plasmas, as for example in the Solar photosphere [10][lTj , or be treated as a small 
perturbation and described within the framework of a perturbation theory [T| IT2T4T9] . 

In this paper we will consider the continuous absorption of EM (electromagnetic) 
radiation in dense partially ionized plasma, with electron density iV e ~ 10 19 cm -3 , 
temperature T m 200007^ and atom density N a m N e . Plasmas with similar parameters 
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are of interest from both the laboratory [2"0"ll2"T] and the astrophysical aspect. Here we 
keep in mind the plasma of the inner layers of the solar atmosphere, as well as of partially 
ionized layers of other stellar atmospheres, for example the atmospheres of DA and DB 
white dwarfs with effective temperatures between lOOOOi^ and 200007^ (see [22142"^] ) . 

Due to the exceptional simplicity of the hydrogen atom, this research is starting 
with the hydrogen case. Under the mentioned conditions the continuous absorption of 
EM radiation in hydrogen plasmas are determined by the following radiative processes: 



where e\ is the energy of a photon with wavelength A, n and I - the principal and the 
orbital quantum number of hydrogen-atom excited states, e k and e^- free electrons with 
energies E = h 2 k 2 /2m and E' = h 2 k' 2 /2m respectively, m the electron mass and H - 
Plank's constant. 

In the considered region of N e and T the processes of electron-if + and electron- 
ic^) inverse bremsstrahlung, photo-detachment of the negative ion H~(ls 2 ), H(ls)- 
H + absorption charge exchange and photo-dissociation of molecular ion i7^(lE+), which 
are described by Eqs. (J2J), (G2) and (TJ1), can be treated in the same way as in previous 
papers [251I26] . Therefore, photo-ionization processes (JT]) will be in the focus of attention 
in the next section, and here let us note only that in this field methods obtained by 
extrapolation (with minor modifications) of the methods developed for weakly and 
moderately non-ideal plasmas [20, [271128] have been used for such processes up until now. 
Moreover, for determination of absorption coefficients characterizing their influence in 
the region of N e > 10 19 cm~ 3 methods based on Cramer's approximation [2T| |2"9] have 
been used so far. So, developing a new modeling way of describing the continual 
absorption in dense partially ionized plasmas, which is the main aim of this paper, 
is still an actual task. 

We take as the landmarks hydrogen plasmas which were experimentally studied 
in [28]: with N e = 6.5 ■ 10 18 cm~ 3 and 1.5 ■ 10 19 cm" 3 , and T = 1.8 ■ 10 4 K and 
2.3 • 10 4 K respectively. The presented modeling way is tested within the optical range 
of photon wavelengths: 300rwTi < A < 500nm. A key parameter of the used numerical 
procedure is determined empirically. It is shown that this procedure already allows 
for determination of the spectral absorption coefficients characterizing all the relevant 
absorption processes in dense partially ionized hydrogen plasmas, at least in the regions 
5 • 10 18 cm" 3 < N e < 5 • 10 19 cm~ 3 and 1.6 • 10 4 K < T < 2.5 • 10 4 K. 

The material of this paper is distributed over Sections 2 and 3. In Sec. 2 the 
following is presented: the approximation of the cut-off Coulomb potential, together 



£ X + H*(n,l) ^H + + e k ,, 
e x + e k + H + -> e k > + H + , 
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with the reasons for applying just this approximation in the considered regions of N e 
and T; the way of obtaining all the partial spectral absorption coefficients, as well as 
their corresponding final expressions. In Sec. 3 the results of the calculations of the 
partial and the total absorption coefficients are presented and discussed. 

2. Theory: the spectral absorption coefficients 

2.1. The approximation of the cut-off Coulomb potential 

Let us note that numerous papers dedicated to development of rather rigorous methods 
of ground-state and exited-atom photo-ionization processes, were published over the 
past few decades (see for example [271I301 - I36] ). Although some of these methods were 
tested on the examples of isolated H-like and He atom systems, the main results of 
those papers are intended for more complex atoms and for atoms in extreme EM fields. 
Therefore we will consider here that an isolated hydrogen atom H*(n,l) should be 
described in the same way as in [37J. The appearance of these papers has certainly 
encouraged intensive investigation of the atom ionization processes in plasmas (see 
for example [38H15] ); however, the regions of temperature and electron density which 
are considered here (N e ~ 10 19 cm -3 and T 20000-/^) have not been examined yet. 
Namely, the results obtained in those papers refer to the regions of extremely large 
electron density and temperature, where the non-relativistic way of describing plasma 
becomes invalid. Besides, the models of ion screening which are used in the mentioned 
papers are not generally appropriate (as it is being shown below) to the considered 
regions of plasma parameters. That is why we examine a new modeling way here, as 
was already mentioned above. 

In the considered hydrogen plasmas processes Eq. (CQ) will present the main 
difficulties in describing the continuous absorption, since under the stated conditions 
the energy of interaction of an excited atom with its neighborhood reaches the order 
of the corresponding ionization potential. It means that these processes have to be 
considered a result of radiative transitions in the whole system "electron-ion pair + the 
neighborhood", namely 



known, many-body processes can sometimes be simplified by a transformation to the 
corresponding single-particle processes in an adequately chosen model potential. 

In this context we will introduce into the consideration the ion Wigner-Seitz radius 
r s -i, given by 



where iVj is the H + ion density, and the characteristic length r c is defined by the relation 





(6) 
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where the quantity U p . c has the meaning of the mean potential in which a free electron 
moves in the considered plasma. Then, we will take into account the fact that, according 
to the results of [16] . the characteristic length r c = r s -i in the hydrogen plasmas which 
are considered here. It means that in these plasmas the radius of the zone where an ion 
and an electron can be treated as an electron-ion pair in Eq. (jSJ) is just r c . 

Here, we will follow the previous papers [17], where it was noted that an adequate 
ion screening potential has to satisfy two main requirements: 

a) it has to be practically equal to Coulomb potential (— e 2 /r) in the mentioned 
zone, i.e. in the region r < r c , 

b) outside of this zone (r > r c ) it has to be practically equal to the mean potential 



where r is the distance from the considered ion, and e - the modulus of the electron 
charge. Condition a) is due by the fact that the averaged potential, created by the rest 
of the considered plasma inside the mentioned zone, where the electron of (H + + e)- 
pair is localized, is equal to zero in the first approximation, since the rest is neutral 
as a whole. As for condition b), it reflects the fact that under such a condition a 
realistic, physically acceptable asymptotics is realized for the electron wave function 
characterizing (H + + e)-pair in the considered plasma, which cannot be Coulomb-like 
as in an isolated electron-ion pair. 

As an adequate model potential for hydrogen plasma with N e ~ 10 19 cm -3 we 
choose, as in [17J, the screening cut-off Coulomb potential, which satisfies conditions a) 
and b), and can be presented in the form 



where the cut-off radius r c is defined by relation (J7J, as it is illustrated by Fig. [T]a. 
This potential was first introduced to the plasma physics in [48] and its properties were 
investigated later in [491 150] . 

In connection with such a choice of the ion screening potential let us note the 
following: the argumentation from [47] is taken into account here, that the often used 
model of Debye-Hiickel (DH) potential is not adequate for description of an electron-ion 
pair in dense non-ideal plasma. It is important that we focus here only on the physical 
meaning of DH potential, leaving out of the consideration all of its disadvantages, which 
have been noted and discussed in [461I5T] ; namely, as it is well known, DH potential has 
the meaning of an average electrostatic potential created by the observed ion and all 
the charged particles from the rest of the plasma. This potential was introduced in [52] 
in order to describe (in accordance with its meaning) some of the thermodynamical 
characteristics of free-particle systems with Coulomb interaction. From this aspect, its 
applications have been justified until now, if the accuracy of the calculations is not very 
important [46] . 

Let us note also that in order to reduce errors which are generated by the application 



U, 




< r < r, 



r c < r < oo 



(8) 
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of DH potential, instead of DH screening radius r D = [kT/{8nN e e 2 )}^ 2 for two- 
component (electron-ion) plasma, another radius r e is often taken, given by 

r e = [fcT/(4 7 riV e e 2 )] 1 /2, (9) 

which in principle characterizes the corresponding single-component system. In |46j it 
was shown that replacing ro with r e is really justified in weakly non-ideal plasmas, 
since the adequate screening radius is approaching just r e when the non-ideality degree 
is approaching zero. Here, it is necessary to remind that by replacing DH screening 
radius rjj with the electron screening radius r e DH potential loses its physical meaning, 
becoming one of the numerous DH-like potentials. 

However, due to its physical meaning, DH potential could not be applied in principle 
to describing electron-ion scattering in dense plasmas. Namely, in such plasmas the 
electron in the considered electron-ion pair is itself a part of the corresponding ion's 
"screening cloud" and thus takes part in the creation of DH potential. This fact is often 
being forgotten, but neglecting it yields an absurd situation in the case of such plasmas 
where DH screening radius ro is approaching the ion Wigner-Seitz radius r s; j, given by 
Eq. El In such a case the complete ion "screening cloud" consists practically of one 
electron from the considered electron- ion pair. 

Just this kind of situation is realized in a plasma of the kind treated in this paper, 
i.e. with N e of the order of magnitude of 10 -19 cm -3 and T of about 20 000K, where the 
ion Wigner-Seitz radius r s; j is close to the ion screening radius r e given by (J9j), making 
an application of DH potential completely unacceptable. Let us emphasize the fact that 
the same is valid for any other potential which has a similar meaning as DH. 

For all these reasons we consider that DH screening potentials, which were used 
in the above mentioned papers [381443] cannot be accepted for an application in any 
region of plasma parameters. As for DH-like potentials, one of them could be accepted 
in principle for an application in some regions of N e and T, but only under the following 
condition: 

-the screening radius in the expression for the used potential guarantees that the 
values of the basic parameters of the considered problem (number and energies of the 
bound states, elastic scattering amplitude for the angle equal to zero as a function of 
the free electron energy, etc.) are very close to the corresponding experimental data. 

This condition could probably be satisfied on the condition of existence of reliable 
experimental data which could be used for determination of the effective screening 
radius by means of the corresponding fitting method, and some rigorous procedure being 
included for describing the bound and free states of the electron in the considered DH-like 
potential. However, for the considered plasmas there are no such reliable experimental 
data in the literature. 

Because of all the things mentioned we consider that the choice of the potential 
U (r; r c ), which is given by Eg. ((8]), is the best solution in the case of the considered dense 
partially ionized plasmas. Besides, an additional advantage of this model potential is 
possibility of obtaining all the needed final results in a compact analytical form. 
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As in |47j, we will take the value U p;c as the zero of energy. After that, the potential 
Eq. ([8]) is transformed to the form 



lUr) { r + r c > 0<r < r c (10) 
0, r c < r < oo, 

which is used in further considerations, and is also illustrated by Fig. [Tib. It is important 
that, through characteristic lengths r S]i and r e defined by Eqs. and ([9]), the cut-off 
radius r c is determined here as a given function of N e and T. Namely, taking that 

r c = a c -r e , (11) 

we can directly determine quantity a c as a function of the ratio r s; j/r e , on the basis of the 
data from [16] about the mean potential energy of an electron in singly ionized plasma. 
We remind that in [16] a quantity is determined which is equal to U p;c /Ud = r D /r c 
(being just designated differently), where Up = —e 2 /r D is the mean potential energy 
determined within DH method [52]. Keeping in mind that r^> = r e /2 1//2 , and taking 
r c /(y/2 ■ ro) = a c , we obtain a c = r c /r e , and then the expression ffTT]) for r c too. As, 
quantity ro/r c is being determined among else in [36], just for the case of singly ionized 
plasma, so the curve in Fig. 2 of this paper is obtained routinely from the curve in Fig 
4 in [16], which represents quantity ro/r c as a function of ratio r s ^/r e . The behaviour 
of a c in a wide region of values of r s; j/r e is presented in FigfJJ 

2.2. Processes (T7]) absorption coefficients. 

For the sake of further considerations it is necessary to have energies E n j of the bound 
states | n, / > which are possible in potential U c (r) and the corresponding radial wave 
functions R ni i(r), as well as the radial wave functions Re,i{t) which correspond to the 
free states \E,l >. These quantities are determined from Schrodinger equation 



d 2 R(r) 2 dR 
— H h 

dr 2 r dr 



*». M(r)) _<«+!> 



h 2 r 



2 



0, (12) 



where W and R(r) denote either E n ^ and R n> i(r) or E and Re,i{t). Here we will find the 
radial wave functions in the form: R n j{r) = P n ,i/ r and Re,i{t) = PE,i/ r , where P n> i{ r ) 
and Pe,i{ t ) can be expressed through the well known analytical functions given in |53j, 
namely: 

- in the region < r < r c all P n ,i(r) and Pe,i(t) with E < e 2 /r c are expressed 
through the corresponding Witteker's functions, and Pe,i{t) with E > e 2 /r c - through 
Coulomb-like ones which are regular at the point r = 0; 

- in the region r c < r < oo all P n ,l{ r ) are expressed through modified Bessel 
functions, and all Pe,i(t) - through spherical Bessel functions. Let us note that the 
bound-state energies (in region W < 0) and the free-state phase shifts (in region W > 0) 
are determined from the condition of continuity of R(r) and dR(r) / dr at the point r = r c . 
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Since the considered photo-ionization processes can be described in the dipole 
approximation (r c << A), the corresponding cross sections, for the non-perturbed bound 
states in potential U c (r), are given by the expressions from [5lj, namely 

a ph (\; n, I, E n>l ) = + S lmax ( J P n,i rP E,i'dr \ , (13) 

where k ph = e x /hc, l max = max(l, I'), and E = E n<l + e x . 

As it is well known, the spectral absorption coefficient K^ h \\; N e , T), characterizing 
the photo-ionization processes in the case of non-perturbed energy levels in the 
potential U c (r), is given by the expression 

K$(\;N e ,T)=^N rh i-a ph (\-,n,l,E n> i)-f st (\,T), n>2, (14) 

n,l 

where a p h(X;n,l, E n j) is the partial photo-ionization cross-section defined by Eq. [131 
N n j is the density of atoms H*(n, I), and factor f st , given by 

f st (X,T) = l-exp(-27rhc/X), (15) 

describes the influence of the stimulated emission. One can see that this expression is 
similar to the one for diluted hydrogen plasma (see for example \1Q\) and, in accordance 
with what was said above, it cannot model the absorption coefficients of dense non-ideal 
plasmas described in [28]. In order to determine the absorption coefficients K p h(X; N e , T) 
which can be applied to the mentioned modeling, we have to take into account additional 
details of the atom-plasma interaction, beside those which are already described by the 
shape of potential ( 1101) . It is known that within the usual way (i.e. as is being done in 
the cases of weakly non-ideal plasmas) this additional influence should be characterized 
by shifts and broadenings of the bound-state energy levels of the considered atoms. 
However, this way is generally applicable in any case, including the case of strongly 
non-ideal plasma too. 

This is confirmed by the approach of [55], main feature of which is treatment of 
electrons in existing atoms H*(n,l) as a Fermi gas of particles, which move in a self- 
consistent field created by immobile ions and other electrons. One of the results of this 
approach is a description of the relevant elementary event in such a gas, i.e. annihilation 
of an electron localized at the j-th proton in the chosen u-th state with a simultaneous 
creation of an electron in some free state, caused by the absorption of a photon with 
energy E\. Namely, it was found that such an event can be described in terms of the 
corresponding probability density, which is practically equal to zero outside of a finite 
interval of free-state energies. This result, despite the fact that it has only qualitative 
significance, suggests that K p h(X; N e ,T) can be obtained by a modification of Eq. (Till) 
based on an adequately chosen probability density p n ,«( e ) °f the perturbed atom energy 
levels with given n and I, characterized by the corresponding shifts A n j and broadenings 
8 n> i. It is assumed that energies e of the perturbed atomic states are dominantly grouped 
around energy e^ ax) = E(n, I) + A„ jZ , inside the interval (e^ ax) - 5 n j/2, e^ ai) + 5 n j/2), 
similarly to the known cases (Gaus, Lorentz, uniform etc.). 
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All this justifies a semi-empirical approach to the considered problem, until a 
corresponding strict method is developed. In this work any (n, /)-shell (with given 
n) of the perturbed hydrogen atom is characterized by only two quantities, namely 
an averaged shift A n and broadening 5 n , which are treated as empirical parameters. 
Consequently, here it is considered that: A n ^ = A n , 5 n j = 5 n and p n ,z(e) describe the 
corresponding uniform distribution of the perturbed energy levels, i.e. 

j_ i (mm) _ i < Sn 
PM={ 5 " "'(max) , 2 L = E n,l + ^n- (16) 

0, lei, -el > — , 
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Here we keep in mind that the uniform distribution can approximate many other 
distributions well (Gaussian, cupola etc.) 

Let us note that although A n is treated here as an empirical parameter, it is 
possible to describe its qualitative behaviour as a function of iV e . Namely, for well-known 
physical reasons all shifts A n j, and consequently A n , have to change proportionally to 
the density of the perturbers, the relative atom-perturber velocity and the characteristic 
perturbation energy. Consequently, we will have it that 

A n = Const. ■ N e ■ v ea (T) ■ e 2 /l(N e , T), (17) 

where v ea (T) and l(N e , T) are the characteristic electron-atom velocity and distance. On 
the basis of the results of [16] in the considered cases (iV e ~ 1 • 10 19 cm -3 , T ~ 2 • 10 4 K) 
any relevant characteristic length has to be close to the radius r e , which is given by 
Eq. ©. Since v ea {T) ~ {k B Tfl 2 and r e ~ (k B T/N e ) 1 / 2 , from ([T7J) the relation follows 

A n = Const. ■ iV e 3/2 , (18) 

which is also in accordance with [55], and will be particularly significant in further 
considerations. 

Here, we will describe the perturbed atomic states in the first order of the 
perturbation theory and, in accordance with what was said above, we will have it that 



K ph (\;N e ,T) 



H ' SnJ e x + e 



fst(X,T), (19) 



Til 

where n > 2, e = S™ x) - 5 n /2, e + = e^ ax) + 8 n /2, e ( ™ x) = E n , t + A n , and 
&ph{^ ■ 6 j n , ^ En,l) is given by Eq. (PC?]) for A^^ = A • E\/(e\ + e), i.e. for the wavelength 
of a photon with energy (e A + e). 

2.3. Processes (H)) ; |3|) and absorption coefficients 

The spectral absorption coefficients characterizing the electron-ion absorption process 
([2]), electron-atom absorption processes ([3]) and ion-atom absorption processes (0J are 
determined here as in the previous papers [25l[26j[56] , dedicated to the same absorption 
processes in some laboratory and astrophysical plasmas. 
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So, the spectral absorption coefficients K ei (X), n ea (X) and Ki a (X) are defined by 
expressions 

K ei {X) = K ei (X,T) -Ne-Nr f st (X,T) S K ei (X,T) ■ N 2 ■ f st (X,T), (20) 
^ea(A) = a pM • S- 1 -N e -N a + K ea (X,T) ■ N e ■ N a ■ f st (X,T), (21) 
A^ e • N a \ A (mkT\ 5 e: 



s -^±^,= 4 fel (22) 

Kia(A) = K ia (X,T) • iV, • iV a = K ia (X, T) -N e -N a - f st (X,T), (23) 

where iV a , and N neg i are the densities of H(ls) atoms, H + ions and negative H~(ls 2 ) 
ions respectively, "eq" denotes that the quantity S ea is determined under the condition 
of thermodynamic equilibrium of the considered system, E~[ is the energy of H~(ls 2 ) 
ion dissociation, and factor f st (X,T) is given by expression (fl5|) . 

The spectral coefficient K e i(X,T), given in [cm 5 ], is calculated by means of the 
corresponding expressions from [53], with the Gaunt factor determined in [57]. The 
spectral coefficient K ea (X,T), given also in [cm 5 ], and the cross-section for photo- 
dissociation of ion H~(ls 2 ) are determined on the basis of the expressions from [58] 
and [59]. Finally, the spectral coefficient K ia (X,T) is determined in the quasi-static 
approximation, which is described in detail in [56], and is taken in the form 

. 12 C(fl A )(fl A /<»o) 4 ( Ui(Rx 



JMA, T) = 0.62 x 10- — — ; ■ «p [-^f- ) , (24) 

where U\(R) is the energy of the electronic ground state 1S+ of molecular ion if^, as a 
function of internuclear distance R, R\ - the resonance internuclear distance for given 
A, ao - the atomic unit of length, and C(R\) - a dimensionless coefficient which is close 
to one. Parameters R\, Ui(R\) and C(R\), as functions of A, are tabulated in [25|l56]. 
The relative contributions of the partial channels, i.e. H(ls)-H + absorption charge 
exchange and molecular ion if^(lS+) photo-dissociation, are obtained by multiplication 
of K ia (X,T) by factors X(z) and [1 — X(z)] respectively. In the considered region of A 
factor X(z) = T(3/2; z)/T(3/2), where z = -U x {R x )/kT. 



3. Results and discussion 



3.1. The characteristics of the cut-off Coulomb potential 

In this paper the approximation of cut-off Coulomb potential ( TlOl) is applied to modeling 
the spectral absorption coefficient obtained in [28] in two experiments with hydrogen 
plasmas, which are treated as a short and a long pulse respectively. In the first case 
(short pulse) plasma with N e = 1.5 ■ 10 19 cm -3 and T = 2.3 ■ 10 4 K was studied, 
while in the second case (long pulse) - it was plasma with N e = 6.5 • 10 18 cm -3 and 
T = 1.8 • lO 4 -^. In the experiments described in [28] plasmas with electron densities up 
to ~ 10 19 cm -3 were created by pulse discharge in quartz capillary. Diagnostics of the 
plasma was carried out on the basis of optical measurements (at A = 632. 8nm), taking 
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into account radial inhomogeneity of the plasma column. The temperature profile is 
defined from independent measurements of brightness and transparency at different 
distances from the center of the capillary. A detailed study is performed just for the 
two above mentioned examples. 

On the basis of Figj2] and Eq. (TTTjl it was found that the cut-off radius r c is equal 
to : 44.964 a.u. for the short pulse, and 55.052 a.u. for the long one. For these values 
of r c the solutions of Eq. f fl2l correspond to the energies of the realized bound state, 
which are presented in Tab's 1 and 2 respectively. The corresponding partial photo- 
ionization cross-sections <r(A; n, I, e n j) are obtained by means of Eq. (fl3|) for n, I and e nj / 
given in Tab's 1 and 2. The behavior of these cross-sections is illustrated in Figj3]and 
FigJUby the examples of photo-ionization cross- sect ions of all realized states with 1=0. 
These figures show qualitative similarity of behavior of the cross-sections in the cases 
r c = 44. 964a. u. and r c = 55. 052a. u. and domination of the cross-sections with n=2. 
One can see a significant difference between the maximal values of the cross-sections 
with n=2 (about 2.70 a.u. and 0.75 a.u.) which correspond to these cases. This fact 
reflects the tendency of a significant decrease of the maximal values of the cross-sections 
for n=2 with an increase of the cut-off radius r c . 

Let us remind that Tab's 1 and 2 characterize the bound states of an electron in 
the potential U c (r) with the values of cut-off radius r c given above. The energies of the 
corresponding ground states approach the value of —Ih, where Ih = 13.598e\^ is the 
tabulated value of the isolated hydrogen-atom ionization potential (see for an example 
NIST Atomic Spectra Database), only when r c — > oo, i.e. when N e — > or T — > oo. 
Also, the energies of the ground states (for the electron densities and the temperatures 
observed) would be close to (—Ih) in the case where instead of U c {r) the potential 
U(r; r c ) would be used (see Fig. 1). 

For each of the considered cut-off radii r c the existence is assumed here of a 
Boltzmann's distribution of the populations N n j of the bound states, given in Tab 
1 or 2, which exist in Eq. (fl9|) . Such a distribution is determined by the corresponding 
values of the total density of neutral hydrogen atoms N a and the temperature T. In 
accordance with [28] here it is taken that: iV a = 1.9 ■ 10 19 cm~ 3 and T = 22980.fr for 
r c = 44.964a.M., and N a = 3.4 • 10 19 cm" 3 and T = 17960^ for r c = 55.052a.w.. As one 
of the consequences, we have it that the total populations of groups of the states with 
same n are equal to: 4.4 • 10 17 , 3.82 • 10 17 and 4.86 • 10 17 cm -3 for n = 2, 3 and 4 in the 
first case, and 1.87 ■ 10 17 , 1.24 ■ 10 17 , 1.44 ■ 10 17 and 1.85 ■ 10 17 cm~ 3 for n = 2, 3 , 4 and 
5 in the second case. 

3.2. The absorption coefficient: the results of the calculations 

In order to compare the obtained theoretical results with the experimental data from 
[28], we have to take into account all the absorption processes which cannot be neglected 
in the considered hydrogen plasmas, i.e. the processes described by Eqs. (JTJ) , ([2]) , (j3j) 
and (jl]). Therefore, when comparing our theoretical results with the experimental data 
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from [28] we have to use the corresponding total spectral absorption coefficient K tat (X), 
namely 

«tot(A) = Kph(X) + Kadd(X), (25) 

where K p h(X) = K p h(X;N e ,T) is given by Eq. ( jl~9l) . and the member n a dd(X) = 
K a dd(X; N e , N a ,T) characterizes the contribution of absorption processes ([2]), ([3]) and 
(jl]). Consequently, we have it that 

Kadd(X) = K ei (X) + n ea (X) + K ia (A), (26) 

where K ei (X) = n ei (X; N e , T), K ea (A) = K ea (A; N e , N a , T) and K ia (A) = K ia (X; N e , N a ,T) 
are the partial spectral absorption coefficients, which are given above by Eqs. fl20l) . fl2T|) . 
([22]), (E3D and (El. 

In accordance with the aims of this work the calculations of the total absorption 
coefficients K to t{X) have been performed for both cases (short and long pulse) in a wide 
region of values of shifts (A n ) and broadenings (5 n ) of atomic levels with n > 2. The 
calculations of Ktot{X) cover the wavelength region 300nm < A < 500nm. However, 
let us emphasize the fact that the values of the experimental absorption coefficient 
n exp (X) characterize not only the bound-free (photo-ionization) processes ([1]) and the said 
additional absorption processes (EJ) - (Hj), but also the bound-bound (photo-excitation) 
processes, which are not considered in this work. Consequently, for the purpose of this 
work the region A < 450nm in the case of short pulse, and A < 425nm in the case of 
long pulse, (see Tabs. 1 and 2) has the real significance, where the considered photo- 
ionization processes dominate in comparison with photo-excitation ones. The results 
of calculations are shown in Figs. 5-10 together with the corresponding experimental 
values Kexp(X) of the spectral absorption coefficient from [28J. 

Figures 5, 6, 7 and 8 illustrate the results of the calculations of K tot (X) in the case 
when A n = const., while Figs. 9 and 10 present the calculations of K iot (A) in the case 
when A n decreases (relative to A 2 ) with increasing n. 

The bottom and the top curves of ntot{X) in Figs. 5 and 6 illustrate strong 
influence of A n on the calculated total absorption coefficient: A n = 5 n = 0.30e^ and 
A„ = 5 n = O.QOeV for the short pulse; A n = 5 n = 0.05eV and A n = 5 n = 0.20eV for 
the long pulse. The groups of three curves, which lie between the corresponding bottom 
and top curves, demonstrate relatively small influence of 5 n on the calculated values of 
K t ot(X): A n = 0A5eV and 5 n = OAOeV, 5 n = 0A5eV and 5 n = 0.50eV for the short 
pulse; A n = 0.125e^ and 5 n = 0.120eV, 5 n = 0.125eV and 5 n = 0.130eK for the long 
pulse. 

Also, the dashed curves on Figs. 5 and 6 demonstrate the behaviour of the spectral 
absorption coefficient K a dd(X), defined by Eqs. ( 1261) and ( 1201 - (1241 . in the case of short 
and long pulses respectively. One can see that the total contribution of electron-ion, 
electron-atom and ion-atom absorption processes, described by Eqs. ([2]), ()3]) and (jl]) 
indeed cannot be neglected in the considered cases. 

Figures 7 and 8 show the curves of n to t(X) calculated with the values of A n and S n 
which are treated as the optimal ones: A n = 0A55eV and 5 n = 0.625eV^ for the short 
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pulse, and A n = 0.13eV and S n = O.lleV for the long one. In order to estimate the 
possible error due to such a choice of shifts and broadenings, the results of calculations 
are shown in Figs. 9 and 10 in the case when A n decreases with the increase of n, 
proportionally to the ionization energies of the corresponding atomic levels. Let us note 
the fact that calculated curves presented in these figures correspond to the optimal values 
of A n=2 and <5 n=2 . One can see that the calculated curves in Figs. 9 and 10 are very close 
to the calculated curves in Figs. 7 and 8, respectively. This fact is reflected in the values 
of A n=2 and <5 n=2 which correspond to the curves in Figs. 9 and 10: A n=2 = 0.49eV^ and 
5 n=2 = 0.65eV for the short pulse, and A„ =2 = 0.14eV and S n= 2 = 0.12eV for the long 
one. 

Beside the curves K fot (A) and K e2 , p (A), the curves K pth (X) are also presented in 
Figs. 7 and 8, which were obtained in [28] for the short and long pulse on the basis of 
the perturbation theory used in [SHIEH]. For the sake of correct interpretation of the 
presented data, let us note the fact that the values of n p .th.W> similarly to K exp (A), 
characterize not only photo-ionization, but also the bound-bound (photo-excitation) 
processes, which are not considered in this work. One can see that in the case of strongly 
non-ideal plasmas (short pulse, N e = 1.5 • 10 19 cm -3 ) the difference between K e:C p(A) and 
Kp.th.W i s so large that it justifies any effort towards development of an alternative 
method of calculation of the strongly non-ideal plasma absorption coefficient. However, 
in the case of long pulse (N e = 0.65 ■ 10 19 cm~ 3 ) the considered plasma is located by its 
parameters in the lower part of the region of strong non-ideality, where the perturbation 
theory should give much better results. This fact is reflected in a significant reduction 
of the difference between K exp (\) and K v th (A). 

Therefore it is important to check whether relation Eq. ( TT8l is valid also for 
N e close to 0.65 • 10 19 cm -3 . Since in the case of constant shifts A n = 0A55eV 
and 0.130eK for short and long pulses respectively, validity of Eq. f fl8|) means that 
0.455/0.130 = (1.5/0.65) 3 / 2 , which is satisfied with an accuracy better than 1%. In the 
case of variable shift we have it that A n=2 = 0.49eV^ and 0.12eV for the short and long 
pulse respectively, and validity of Eq. ({TBI means now that 0.49/0.14 = (1.5/0.65) 3//2 , 
which is satisfied with the same accuracy. 

This fact offers a possibility to determine A n or A n=2 for any N e and T from 
intervals 0.65 • 10 19 cm" 3 < N e < 1.5 • 10 19 cm- 3 and 1.8 • 10 4 K < T < 2.3 • 10 4 K, and 
probably in some significantly wider regions, at least for 5-10 18 cm -3 < N e < 5-10 19 cm~ 3 
and 1.6- 10 4 -fT < T < 2.5- 10 4 i^. Then, since it has been established that the influence of 
S n is significantly weaker than the influence of A n , we can determine 5 n for any N e from 
those intervals, taking 5 n = A n , as in the examples illustrated by the dashed curves on 
Figs. [7] and E 

On the grounds of all that has been said one can conclude that the presented 
method can already be used for calculations of the spectral absorption coefficients 
of dense hydrogen plasmas in the regions N e ~ 10 19 cm -3 and T e pa 2 • 10 4 i^, as 
long as electron-if + and electron- H(ls) inverse bremsstrahlung, negative ion H~(ls 2 ) 
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photo-detachment, H(ls)-H + absorption charge exchange and molecular ion i?2~(lE+) 
photo-dissociation can be described as in this paper. Let us note the fact that with 
some modifications related to the atom photo-ionization processes, which should enable 
description of the influence of the atom core, the presented method can be applied 
to other kinds of laboratory dense hydrogen-like plasmas. Most of all, we mean 
alkali metals, helium and other rare-gas plasmas. Also, the obtained results can be 
of interest for some astrophysical plasmas, namely the plasma of the inner layers of 
solar atmosphere, as well as the plasmas of partially ionized layers of some other stellar 
atmospheres (for example some DA and DB white dwarfs). 

Further development of the method described requires first of all an improvement 
of the procedure used, in order to replace the semi-empirical parameters with ones 
determined within the procedure itself. Another step would be including into 
consideration atom bound-bound (photo-excitation) processes which have been omitted 
here, and extending the region of the method's applicability to the long wavelengths. 
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Figure 1. Potentials U(r;r c ) and U c (r), where r c is the cut-off parameter. 



Table 1. The energies of the bound states in the potential U c (r) in the case of short 
pulse (cut-off radius r c = 44.964 a.u.): the principal and orbital quantum numbers n 
and I, and the corresponding energies E n j in [cm -1 ]. 
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Figure 2. The parameter a c = r c /r e as the function of the ratio r s -i/r e , where r s -i 
and r e are given by Eqs. © and ©. 



Table 2. The energies of the bound states in the potential U c (r) in the case of long 
pulse (cut-off radius r c — 55.052 a.u.): the principal and orbital quantum numbers n 
and and the corresponding energies E 7li i in [cm -1 ]. 
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Figure 3. The photo-ionization cross-section for the bound state with I = and 
n = 2, 3 and 4 in the potential U c {r) for short pulse (r c = 44. 964a. u.). 
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Figure 4. The photo-ionization cross-section for the bound state with I = and 
n = 2, 3, 4 and 5 in the potential U c (r) for the long pulse (r c = 55. 052a. u.). 
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Figure 5. The influence of the shift and broadening to the calculated total 
absorption coefficient. Short pulse: 1 - absorption coefficient Kt t(A) calculated in 
the approximation of the constant shift with A„ = 5 n = 0.3eV; 2- the same, but with 
A„ = S n = 0.6eV; 3, 4 and 5 - the same, but with A„ = 0A5eV and S n — OAeV, 
0A5eV and 0.5eV respectively. Dashed curve - absorption coefficient K a dd(A) which 
characterizes the total contribution of absorption processes © - ((U). Bold curve - the 
experimental absorption coefficient from |28j . 
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Figure 6. The influence of the shift and broadening to the calculated total 
absorption coefficient. Long pulse: 1 - absorption coefficient K tot (A) calculated in the 
approximation of the constant shift with A„ = S n = 0.05eV^; 2 - the same, but with 
A„ = S n = 0.2eV; 3, 4 and 5 - the same, but with A„ = 0.125eF and S n = O.UeV, 
0.125eV and 0.13ey respectively. Dashed curve - absorption coefficient K a dd(A) which 
characterizes the total contribution of absorption processes © - ((U). Bold curve - the 
experimental absorption coefficient from |28j . 
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Figure 7. Short pulse: 1 - absorption coefficient K to t(X) calculated in the 
approximation of the constant shift with the optimal values of A„ and 5 n , namely 
A„ = 0.455eV and S n = 0.525e^; 2 - absorption coefficient K{ G t(A) calculated in the 
same approximation with A„ = 5 n = 0A55eV. Bold and doted curve - the absorption 
coefficients obtained in j28) experimentally and by means of the perturbation theory. 
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Figure 8. Long pulse: 1 - absorption coefficient K tot (X) calculated in the 
approximation of the constant shift with the optimal values of A„ and S n , namely 
A„ = O.lSeV and 5 n — O.lleT^; 2 - absorption coefficient Kf t(A) calculated in the 
same approximation with A„ = 5 n = 0.13eV^. Bold and dashed line - the absorption 
coefficients obtained in j28) experimentally and by means of the perturbation theory. 
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Figure 9. Short pulse: thin line - absorption coefficient K tot (X) calculated in the 
approximation of the variable shift with the optimal values of A n= 2 and 5 n= 2, namely 
A rl= 2 = 0.49eV" and S n= 2 — 0.55eV^. Bold line - the experimental absorption coefficient 
from gS]. 
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Figure 10. Long pulse: thin line - absorption coefficient K tot (X) calculated in the 
approximation of the variable shift with the optimal values of A n= 2 and S n =2, namely 
A n= 2 = 0.14eV and S n= 2 = 0.12eV. Bold line - the experimental absorption coefficient 
from gS]. 



